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STABILITY OF DENSITIES FOR PERTURBED DIFFUSIONS AND MARKOV 

CHAINS* * ** ” 

V. Konakov\ a. Kozhina^ and S. Menozzi^ 


Abstract. We study the sensitivity of the densities of non degenerate diffusion processes and related 
Markov Chains with respect to a perturbation of the coefficients. Natural applications of these results 
appear in models with misspecified coefficients or for the investigation of the weak error of the Euler 
scheme with irregular coefficients. 

Resume. Nous etudions la sensibilite des densites de processus de diffusion non degeneres et des 
Chaines de Markov associees par rapport a une perturbation des coefficients. Ces resultats trouvent 
des applications naturelles dans I’etude de modeles avec incertitude sur les coefficients ou pour I’analyse 
de I’erreur faible du schema d’Euler a coefficients irreguliers. 
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1. Introduction 


1.1. Setting. 

For a fixed given deterministic final horizon T > 0, let us consider the following multidimensional SDE: 

dXt = b{t,Xt)dt + ait,Xt)dWt, tG[0,r], (1.1) 

where b : [0,T] x —>■ a : [0,T] x R^^ —>■ R'^ ® R'^ are bounded coefficients that are measurable in time 

and Holder continuous in space (this last condition will be possibly relaxed for the drift term b) and IF is a 
Brownian motion on some filtered probability space (H, (J^t)(>o,P). Also, a{t,x) := aa*{t,x) is assumed to 

be uniformly elliptic. In particular those assumptions guarantee that (HID admits a unique weak solution, see 
e.g. Bass and Perkins |BP09) . |Menll) from which the uniqueness to the martingale problem for the associated 
generator can be derived under the current assumptions. 
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We now introduce, for a given parameter e > 0, a perturbed version of dm with dynamics: 

= b,it,X^^^)dt + aS,X^^^)dWt, t e [0,T], (1.2) 

where ; [0, T] x K.'^ —>■ : [0, T] x —>■ ® satisfy at least the same assumptions as b, a and are in 

some sense meant to be close to 5, cr when e is small. 

It is known that, under the previous assumptions, the density of the processes (W)t>o, exists in 

positive time and satisfies some Gaussian bounds, see e.g Aronson |Aro59| . Sheu [She91) or [DM 10) for extensions 
to some degenerate cases. 

The goal of this work is to investigate how the closeness of {bg,ae) and (6, cr) is reflected on the respective 
densities of the associated processes. Important applications can for instance be found in mathematical finance. 
If the dynamics of HU models the evolution of the (log-)price of a financial asset, it is often very useful to 
know how a perturbation of the volatility a impacts the density, and therefore the associated option prices (see 
also Corielli et al. jCFPIOj or Benhamou et al. [BGMlOj for related problems). 

In the framework of parameter estimation it can be useful, having at hand estimators (b^jCrg) of the true 
parameters (b,a) and some controls for the differences \b — 5e|, |cr — ae\ in a suitable sense, to quantify the 
difference — p of the densities corresponding respectively to the dynamics with the estimated parameters and 
the one of the model. 

Another important application includes the case of mollification by spatial convolution. This specific kind of 
perturbation is useful to investigate the error between the densities of a non-degenerate diffusion of type (EH) 
with Holder coefficients (or with piecewise smooth bounded drift) and its Euler scheme. In this framework, 
some explicit convergence results can be found in [KM 15) . 

More generally, this situation can appear in every applicative field for which the diffusion coefficient might 
be misspecified. 

The previously mentioned Gaussian bounds on the density are derived through the so-called parametrix 
expansion which will be the crux of our approach. Roughly speaking, it consists in approximating the process 
by a proxy which has a known density, here a Gaussian one, and then in investigating the difference through 
the Kolmogorov equations. Various approaches to the parametrix expansion exist, see e.g. Il’in et al. [IK062| . 
Friedman [Fri64] and McKean and Singer [MS67j . The latter approach will be the one used in this work since 
it appears to be the most adapted to handle coefficients with no a priori smoothness in time and can also be 
directly extended to the discrete case for Markov chain approximations of equations EH) and EH- Let us 
mention in this setting the works of Konakov et al., see [KMOPj . [KM02] . 

Our stability results will also apply to two Markov chains with respective dynamics: 

+b{tk,Yt^)h +a{tk,Yt,^)Vh^k+i,YQ = X, 

+ b ,{ tk , Yl^^)h + aSk , = X , (1-3) 

where /i > 0 is a given time step, for which we denote for all k > 0, tk '.= kh and the (^fc)fc>i are centered 
i.i.d. random variables satisfying some integrability conditions. Again, the key tool will be the parametrix 
representation for the densities of the chains and the Gaussian local limit theorem. 


1.2. Assumptions and Main Results. 

Let us introduce the following assumptions. Below, the parameter e > 0 is fixed and the constants appearing 
in the assumptions do not depend on e. 
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(Al) (Boundedness of the coefficients). The components of the vector-valued functions b{t,x),be{t,x) and 
the matrix-functions a{t,x),ae{t,x) are bounded measurable. Specifically, there exist constants Ki,K 2 > 0 s.t. 

sup \b{t,x)\ + sup \be{t,x)\ < Ki, 

(i,£c)G[0,T] xR'* (t,£c)G[0,T] xR*^ 

sup |cr(t,ai)|-I- sup \< Je { t , x )\< K2 . 

(t,a:)G[0,T]xR‘^ (i,a:)G[0,T] xR"^ 

(A2) (Uniform Ellipticity). The matrices a := (ja*,as := (JeU* are uniformly elliptic, i.e. there exists 
A>1, V(t,cr,Oe [ 0 ,T]x (R'')2, 

A-'ICP < («(t,a:)^,0 < A|CP, A-i|ep < (ae(t,x)C,e) < Aj^p. 

(A3) (Holder continuity in space). For some 7 G (0,1] , k < 00 , for all t G [0,T], 

|cr(t,x) - a(t,y)\ + \(Te{t,x) - cre(t,y)| < k|x - . 

Observe that the last condition also readily gives, thanks to the boundedness of cr, CTe that a, Oe are also uni¬ 
formly y-Holder continuous. 

For a given e > 0, we say that assumption (A) holds when conditions (A1)-(A3) are in force. Let us now 
introduce, under (A), the quantities that will bound the difference of the densities in our main results below. 
Set for £ > 0: 

AE,b,oo := sup | 6 (t,x) - 6 E(i,a;)l , Vg G (1,-foo), := sup ||5(t,.) - 6 E(i, •)IIl-!(rO- 

(t,a;)G[0,T]xR‘* *6(0,T] 

Since cr, cte are both y-Holder continuous, see (A3) we also define 

Ae,<t ,7 := sup |cr(u,.) - crE(M, .)|^, (1.4) 

iiG[0,T] 

where for 7 G (0,1], |.|-y stands for the usual Holder norm in space on 0 M'^) (space of Holder 

continuous bounded functions, see e.g. Krylov |Kry96|) i.e. : 


I/I7 := sup |/(x)| -f [f]j, [f]^ := 


sup 


\f{x) - f{y)\ 

\x-y\'r 


The previous control in particular implies for all {u,x,y) G [0,T] x 

|a(u, x) - a(u, y) - ae{u, x) + Oe(u, y)\ < 2(^2 -I- K)AE,a. 7 |x - 
We eventually set for g G (1, -foo], 

AE,7,q' .= AE,cr,7 ^e,b^q-; (^A) 

which will be the key quantity governing the error in our results. 

We will denote, from now on, by (7 a constant depending on the parameters appearing in (A) and T. We 
reserve the notation c for constants that only depend on (A) but not on T. The valnes of C, c may change from 
line to line and do not depend on the considered parameter £. Also, for given integers i,j G N s.t. i < j, 
we will denote by [i, j] the set + - ■ ■ , j}- 

We are now in position to state onr main results. 
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Theorem 1 (Stability Control for diffusions). Fix £ > 0 and a final deterministic time horizon T > 0. Under 

(A) and for q > d, there exist constants C := C{q) > 1 and c G (0,1] s.t. for all 0 < s < t < T, {x,y) G 

Pc{t - s,y - x)~^\{p-Pe){s,t,x,y)\ < (1.6) 

wherep{s,t,x, .),Pe{s,t,x,.) respectively stand for the transition densities at timet of equations p.ll) . (11.21) start¬ 
ing from X at time s. Also, we denote for a given c > 0 and for all u > 0, z G , pdu, z) := 

Remark 1 (About the constants). ITe mention that the constant C := C{q) in (11.611 explodes when q f d and 

is decreasing in q (see Lemmas\^ and\M below for details). In particular, it can be chosen uniformly as soon as 
q>qo>d. 

Before stating our results for Markov Chains we introduce two kinds of innovations in (11.31) . Namely: 

(Ig) The i.i.d. random variables {fk)k>i are Gaussian, with law J\f{0,ld). In that case the dynamics in (11.31) 
corresponds to the Euler discretization of equations HU and (HU). 

(Ip. m') For a given integer AI > ‘2d + 5 + y, the innovations {fk)k>i are centered and have density f^ which 
has, together with its derivatives up to order 5, at least polynomial decay of order M. Namely, for all 2 G 
and multi-index u, |u| < 5: 

\D’'fdz)\<CQM{z), (1.7) 

where we denote for all r > d, z G Qr{z) '■= Cr (i+ | ^ | )r ; J^d dzQr{z) = 1. 

Theorem 2 (Stability Control for Markov Chains). Fix £ > 0 and a final deterministic time horizon T > 0. 
For h = T/N, G N* := N \ {0}, we set for i G N, ti := ih. Under (A), assuming that either (Ig) or (Ip^m) 
holds, and for q > d there exist C := C{q) > 1, c G (0, 1] s.t. for all 0 < ti < tj < T, (x, y) G (R.'^)^; 

Xc{t] - ti,y - x)~^\{p^ - p’l){ti,tj,x,y)\ < CAs,j,q, ( 1 . 8 ) 

where p^{ti,tj,x, ■) respectively stand for the transition densities at time tj of the Markov Chains 

Y and in HU starting from x at time ti. Also: 

- If (Ig) holds: 

xSj -ti,y- x) := pdtj -ti,y- x), 

with Pc as in Theorem [TJ 

- Ifilpjci) holds: 

Xc(tj -ti,y- x) := (^_^^)rf/ 2 QM-(d+ 5 + 7 ) [jj-ZrfFpTJf.) ' 

The behavior of C := C{q) is similar to what is described in Remarkl^ 

1.3. On Some Related Applications. 

1.3.1. Model Sensitivity for Option Prices. 

Assume for instance that the (log)-price of a financial asset is given by the dynamics in ( 11 . 11 ) . Under suitable 
assumptions the price of an option on that asset writes at time t and when Xt = x as ]E[/(exp(A.^’^))] up to 
an additional discounting factor. In the previous expression / is the pay-off function. For a rather large class 
of pay-offs, say measurable functions with polynomial growth, including irregular ones, Theorem [Tj allows to 
specifically quantify how a perturbation of the coefficients impacts the option prices. Precisely for a given £ > 0, 
under (A): 



\£dUT,xJ)\ := |E[/(exp(A^"))] -E[/(exp(A^")’‘’"))]| < 


f{exp{y))pc{T - t, x, y)dy. 
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This previous control can be as well exploited to investigate perturbations of a model which provides some 
closed formulas, e.g. a perturbation of the Black and Scholes model that would include a stochastic volatility 
taking for instance (Je{x) = cr + eV'(x) for some bounded y-Holder continuous function '0 and e small enough. 
In that case, assuming that the drift is known and unperturbed, we have oo = We ~ W-y = WWW- 

In connection with this application, we can quote the work of Corielli et al. [CFPlOj who give estimates on 
option prices through parametrix expansions truncating the series. Some of their results, see e.g. their Theorem 
3.1, can be related to a perturbation analysis since they obtain an approximation of an option price for a 
local volatility model in terms of the Black-Scholes price and a correction term corresponding to the first order 
term in the parametrix series. A more probabilistic approach to similar problems can be found in Benhamou 
et al. [BGMlOj . However, none of the indicated works indeed deals with the global perturbation analysis we 
perform here. 


1.3.2. Weak Error Analysis 


It is well known that if the coefficients 6, a in (HD are smooth and a satishes the non-degeneracy condition 
(A2), then the weak error on the densities for the approximation by the Euler scheme is well controlled. 
Precisely, for a given time step h > 0, let us set for i G N,ti := ih. Introduce now, for a given starting time 
ti, the Euler scheme Xf. = x, Vj > i, = X^. + b(tj,X^W^ + a{tj, X^WO^tj+i — W^tj) and denote by 

.) its density at time tj. The dynamics of the Euler scheme clearly enters the scheme (II.3p . It 
can be derived from Konakov and Mammen [KM02j (see also Bally and Talay |BT96j for an extension to the 
hypoelliptic setting) that: 


\p - < Chpc{tj -U,y- x). 

If the coefficients in dHD are not smooth, it is then possible to use a mollihcation procedure, taking for x G 
bs{t,x) := b{t,.) -k pe{x),ae{t,x) := a(t, .) * Pe{x) with pe(x) := e~‘^p{x(e) and p G K+), p(x)dx = 

l,supp(p) C K for some compact set K of For the mollifying kernel pe, one then easily checks that for 
7-H61der continuous in space coefficients b, a there exists C s.t. 


sup \b{t,.) - bs{t,.)\^ < Ce'^, sup \a{t,.) - ■)\r, < Ce'^ '', 776 ( 0 , 7 ). 

tG[0,T] te[0,T] 


(1.9) 


The important aspect is that we lose a bit with respect to the sup norm when investigating the Holder norm. 
We then have by Theorems [T] and [2] and their proof, that, for y-Holder continuous in space coefficients b, a and 
taking q = 00 , there exist c, C s.t. for all 0 < s < t < T, 0 < ti < tj < T, {x, y) G (K'^)^: 

\{P-Pe){s,t,x,y)\ < CCr^£'^~'^Pc{t- s,y - x), Hp’^ -p^){U,tj,x,y)\ < CCr,£'^~'^Pc{tj - U,y-x), 


where the constant explodes when 77 tends to 0 . 

To investigate the global weak error (p — pWiU, tj,x, y) = {{p — Pe) + {Pe — P^) + {Pe ~ tj^x, y), it 

therefore remains to analyze the contribution {p^ —p^){ti, tj,x, y). The results of |KM02j indeed apply but yield 
\{pe — p^){ti,tj,x,y)\ < Cehpc{tj —ti,y — x) where Cg is explosive when e goes to zero. The global error thus 
writes: 

\ip-p’'){U,tj,x,y)\ < C{Cne'^~'^ + Ceh}pc(tj -ti,y-x), 
and a balance is needed to derive a global error bound. This is precisely the analysis we perform in |KM15j . In 
this work, we extend to densities (up to a slowly growing factor) the results previously obtained by Mikulevicius 
and Platen [MP91| on the weak error, i.e. they showed |E[/(At) — /(A.^)] | < provided / G R). 

Precisely, we obtain through a suitable analysis of the constants which respectively depend on behavior of 

the parametrix series and of the derivatives of the heat kernel with mollified coefficients, that \p—p^\{ti,tj,x, y) < 

function ip{h) going to 0 as h —>■ 0 (which is induced by the previous loss of rj 
in (11.91) 1. In the quoted work, we also obtain some error bounds for piecewise smooth drifts having a countable 
set of discontinuities. This part explicitly requires the stability result of Theorems [U Theorems [5] for q < -|-c». 
The idea being that the difference between the piece-wise smooth drift and its smooth approximation (actually 
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the mollification procedure is only required around the points of discontinuity), is well controlled in norm, 
q < + 00 . 

1.3.3. Extension to some Kinetic Models 

The results of Theorems [T] and [5] should extend without additional difficulties to the case of degenerate 
diffusions of the form: 


dX} = b{t,Xt)dt + a{t,Xt)dWt, 
dXf = X^dt, 


denoting Xt = (X^^X^), under the same previous assumptions on b,a when we consider perturbations of the 
non-degenerate components, i.e. for a given e > 0, xj^^ = (xl'^^\ where: 


= b,(t, Xi^^)dt + aS, X^^^)dWu 

dXf’^'"'’ = Xl’^’^Ut. 


( 1 . 11 ) 


Indeed, under (A), the required parametrix expansions of the densities associated with the solutions of equation 
(II.IOL (11.111) have been established in |KMM10) . 


1.3.4. A posteriori Controls in Parameter Estimation. 

Let us consider to illustrate this application a parametrized family of diffusions of the form: 


dXt = b{t, Xt)dt + a{d, t, Xt)dWt, (1.12) 

where S 0 C the coefficients &, a are smooth, bounded and the non-degeneracy condition (A2) holds. 
A natural practical problem consists in estimating the true parameter d from an observed discrete sample 
(Atn)jgjg_„| where the {(ii)ig|o,ra]} form a partition of the observation interval, i.e. if T = 1, 0 = fg < < • < 

f" = 1 ' 

''ra 

Introducing the contrast: 


1 " 

[/"(!?) := -^ [log (det(a(d,fr-i,A*._J)) + {a-\d,tU,Xtj_JXr,Xr) 

77 . L 
i—l 


V*e [I,n],Ar := 


Xt^ - Xt2_^ 

v't? - t?-i 


and denoting by the corresponding minimizer, it was shown by Genon-Catalot and Jacod |GCJ93] that 
under P*^, ^/n{dn — d) converges in law towards a mixed normal variable S which is, conditionally to Pi := 
[ 0 . 1 ]}) centered and Gaussian. For a precise expression of the covariance which depends on the whole 
path of (At)tg[o_i] we refer to Theorem 3 and its proof in [GGJ93] . 

This means that, when n is large, conditionally to Pi, we have on a subset Hen which has high probability, 
that — i?! < -^ for a certain threshold C. Setting ae^{t,x) := cr(i?„,t,x) and with a slight 

abuse of notation a{t,x) := a{'d,t,x), one gets that, on Cl: 

\<j{t,x) - crs^{t,x) - {a{t,y) - ae^{t,y))\ < |a; - ?/| A ^ \<J - cre\r, < (Cn"^/^)^"’', g G (0,1]. 

We can then invoke our Theorem [1] to compare the densities of the diffusions with the estimated parameter and 
the exact one in (I1.12p . 

The paper is organized as follows. We recall in Section [5] some basic facts about parametrix expansions for 
the densities of diffusions and Markov Ghains. We then detail in Section [3] how to perform a stability analysis 
of the parametrix expansions in order to derive the results of Theorems [T] and 
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2. Derivation of formal series expansion for densities 


2.1. Parametrix Representation of the Density for Diffusions 

In the following, for given {s,x) G R’*' x R'^, we use the standard Markov notation {X^’^)t>s to denote the 
solution of (ffH) starting from x at time s. 

Assume that {X^’^)t>s has for all t > s a smooth density p{s, t, x ,.) (which is the case if additionally to (A) 
the coefficients are smooth see e.g. Friedman [Fri64] 1. We would like to estimate this density at a given point 
y G R'^. To this end, we introduce the following Gaussian inhomogeneous process with spatial variable frozen 
at y. For all {s,x) G [0,T] x R"^, t > s we set: 

Xf = a: + f a(u, y)dWu- 

J S 

Its density readily satisfies the Kolmogorov backward equation: 


duP'^iu, t, z, y) + Lypy{u, t,z,y) = 0, s <u <t,z G R”^, 

.,y) 6y{.), 


where for all (p G C'q(R‘^,R) (twice continuously differentiable functions with compact support) and z G R'^: 

= ^Tr (aa*(u,y)Dlip(z)) , 


stands for the generator of X'^ at time u. 

On the other hand, since we have assumed the density of X to be smooth, it must satisfy the Kolmogorov 
forward equation (see e.g. Dynkin |Dyn65|). For a given starting point x G at time s. 


dup(s, u, X, z) — L’^p{s, u,x, z) = 0, s < u < t, z G R'^, 
p{s,U,X,.) -> Sa :{.), 

U^s 


( 2 . 2 ) 


where L* stands for the formal adjoint (which is again well defined if the coefficients in (11.11) are smooth) of 
the generator of (II.ip which for all ip G Cq(R'^,R), z G R'^ writes: 

Lup{z) = ^Tr {aa* {u, z)DIp{z)) + {b{u, z), D^p{z)). 

Equations (12.ip , (12.2p yield the formal expansion below which is initially due to McKean and Singer [MS67] . 

{p-py){s,t,x,y)= duduj dzp{s,u,x,z)py{u,t,z,y) 

Js dRd 

= du dz{dup{s,u,x,z)p^{u,t,z,y)+pis,u,x,z)duP^iu,t,z,y)) 

Js 

= du dz {l*uP{s, u, X, z)p^(m, t, z, y) - p{s, u, x, z)Llp'^{u, t, z, y)^ 


= du dzp{s,u,x,z){Lu-Ll)py{u,t,z,y), 
Js Jr<^ 


(2.3) 


using the Dirac convergence for the first equality, equations (lO) and dUD for the third one. We eventually take 
the adjoint for the last equality. Note carefully that the differentiation under the integral is also here formal 
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since we would need to justify that it can actually be performed using some growth properties of the density 
and its derivatives which we a priori do not know. 

Let us now introduce the notation 


f ® g{s,t,x,y) = du dzf(s,u,x,z)g{u,t,z,y) 
Js JR-i 


for the time-space convolution and let us define p{s,t,x,y) := p'^{s,t,x,y), that is in p(s,t,x,y) we consider 
the density of the frozen process at the final point and observe it at that specific point. We now introduce the 
parametrix kernel: 

H{s,t,x,y) := {L^ - Ls)p{sfi,x,y) := {L^ - LDp^{sfi,x,y). (2.4) 

With those notations equation (12.31) rewrites: 

{p-p){s,t,x,y) =p(S) H{s,t,x,y). 

From this expression, the idea then consists in iterating this procedure for p{s,u,x, z) in ()2.3I) introducing the 
density of a process with frozen characteristics in z which is here the integration variable. This yields to iterated 
convolutions of the kernel and leads to the formal expansion: 


p{s,t,x,y) = '^p(^ H^'"\s,t,x,y), 


(2.5) 


r—0 


where p (gi = p, = H ® r > 1. Obtaining estimates on p from the formal expression (ESI) 

requires to have good controls on the right-hand side. The remarkable property of this formal expansion is now 
that the right-hand-side of (12.51) only involves controls on Gaussian densities which in particular will provide, 
associated with our assumption (A) a smoothing in time property for the kernel H. 

Proposition 1. Under the sole assumption (A), fort > s, the density of solving (11.11) exists and can be 
written as in 

Proof. The proof can already be derived from a sensitivity argument. We first introduce two parametrix series 
of the form (12.51) . Namely, 


p{s,t,x,y) :=p{s,t,x,y) + '^p H^'"\s,t,x,y) 


r—1 


and 


Pe{s,t,x,y) ■.= Pe{s,t,x,y) + ^Pe ® H'f\s,t,x,y). 


( 2 . 6 ) 


(2.7) 


r—1 


Let us point out that, at this stage, p and Pe are defined as sums of series. The purpose is then to identify those 
sums with the densities of the processes at point y. 

The convergence of the series (12.61) and (12.71) is in some sense standard (see e.g. |Menll) or Friedman [Fri64] 1 
under (A). We recall for the sake of completeness the key steps for (12.6|) . 

From direct computations, there exist ci > l,c S (0,1] s.t. for all T > 0 and all multi-index a, |a| < 8, 


VO < w < t < T, (z, j/) e (M ) , \Dfp{u, t, z,y)\< 


Cl 


(t-n)l“l/2 


Pc{t -u,y- z), 


( 2 . 8 ) 


where 


Pcit -u,y- z) = 


^d/2 


(27r(< — n))^/2 


exp 


\y- 


2 t — u 
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stands for the usual Gaussian density in with 0 mean and covariance {t — u)c ^Id- From (12.81) . the bound¬ 
edness of the drift and the Holder continuity in space of the diffusion matrix we readily get that there exists 
Cl > 1, ce (0,1], 

Cl (1 V 

\H[u,t,z,y)\< Pc{t-u,z -y). (2.9) 

Now the key point is that the control (12.9p yields an integrable singularity giving a smoothing effect in time 
once integrated in space in the time-space convolutions appearing in (12.61) and p.7l) . It follows by induction 
that: 


\p® H^'^\s,t,x,y)\ < ((1 + {i - - s,y - x){t - sf 


((1 [r(^)]' 

fO+Tf) 


-Pc{t - s,y - x){t - s) ^ , 


( 2 . 10 ) 


where for a,b > 0,B{a,b) = fgt ^^“(1 — t) stands for the /3-function, and using as well the identity 

B{a,b) = for the last equality. These bounds readily yield the convergence of the series as well as a 

Gaussian upper-bound. Namely 


p{s,t,x,y) < ciexp((l V '^)/^)ci[(t - s)^/^])pc(i - s,y - x). 


( 2 . 11 ) 


An important application of the stability of the perturbation consists in considering coefficients b^ := b * 
Ce,cre := cr*Ce in dm, where Q is a mollifier in time and space. For mollified coefficients which satisfy the 
non-degeneracy assumption (A2), the existence and smoothness of the density of is a standard property 
(see e.g. Friedman |Fri64) . |Fri75) b It is also well known, see |KM00] . that this density p^ admits the series 
representation given in dm- Observe now carefully that the previous Gaussian bounds also hold for p^ uniformly 
in s and independently of the mollifying procedure. This therefore gives that 

Pe{s, t, X, y) —> p(s, t, X, y), (2.12) 

£->■0 

boundedly and uniformly. Thus, for every continuous bounded function / we derive from the bounded conver¬ 
gence theorem and (12.111) that for all 0 < s < t < T, x S 

]gjy(^G).s,^)] ^ / f(jj)p^i^s,t,x,y)dy — f{y)p{s,t,x,y)dy. (2.13) 

Taking / = 1 gives that p(s, t, x, y)dy = 1 and the uniform convergence in (12.121) gives that p(s, t, x ,.) is non 
negative. We therefore derive that p{s,t,x,-) is a probability density on 

On the other hand, under (A), we can derive from Theorem 11.3.4 of [SV79] that (Af ) 5 gfo t] (As);ig[o Tl ■ 
This gives that for any bounded continuous function /: 

E[/(aG).«.-)] E[/(A*’")]. 

6^0 

This convergence and (12.131) then yield that the random variable admits p{s,t,x, ■) as density. 

We can thus now conclude that the processes A, A*^^) in (fTTl) . (HH) have transition densities given by the 
sum of the series ( 1 ^ . dm- 

□ 
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2.2. Parametrix for Markov Chains 

One of the main advantages of the formal expansion in (12.Sp is that it has a direct discrete counterpart in 
the Markov chain setting. Indeed, denote by the Markov chain with dynamics (El starting from x 

at time U. Observe first that if the innovations {^k)k>i have a density then so does the chain at time tk- 

Let us now introduce its generator at time ti, i.e. for all if G C'q(R‘^,R), x G 

L^fix) := h-^E[f{Y,^;f^)-f{x)]. 

In order to give a representation of the density of p^{ti,tj,x,y) of at point y for j > i, we introduce 
similarly to the continuous case, the Markov chain (or inhomogeneous random walk) with coefficients frozen in 
space at y. For given {ti,x) G [0,T] x R^, tj > ti we set: 


.= a- h^/^'Y^a{tk,y)ik+i, 

k—i 

and denote its density by {ti,tj,x ,.). Its generator at time ti writes for all f G C'q(R‘^,R),x G 

L’l:yf{x) = h-^nf{Y,X''^)-f{x)]. 

Using the notation p^{ti,tj,x,y) := p^’'^{ti,tj,x,y), we introduce now for 0 < i < j < fV the parametrix 
kernel: 

H^(ti,tj,x,y) := (LJ* - + h, ,x, y). 

Analogously to Lemma 3.6 in |KM00) . which follows from a direct algebraic manipulation, we derive the following 
representation for the density which can be viewed as the Markov chain analogue of Proposition [TJ 

Proposition 2 (Parametrix Expansion for the Markov Chain). Assume (A) is in force. Then, for 0 < U < 
tj < T, 

j-i 

p^{ti,tj,x,y) = ^P^ H^’^'^Hu,tj,x,y), 

r—0 

where the discrete time convolution type operator i^ih is defined by 

j-i-i j, 

f ®hg{ti,tj,x,y) = '^h f{ti,ti+k,x,z)g{ti+k,tj,z,y)dz. 
k=o 

Also g = g and for all r > := i^k denotes the r-fold discrete convolution of the 

kernel H^. 


3. Stability of Parametrix Series. 

We will now investigate more specifically the sensitivity of the density w.r.t. the coefficients through the 
difference of the series. For a given fixed parameter e, under (A) the densities p(s,t,x, .),Pe(s,t,x, •) at time 
t of the processes in dni), dOj) starting from x at time s both admit a parametrix expansion of the previous 
type. 
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3.1. Stability for Diffusions: Proof of Theorem [T] 

Let us consider the difference between two parametrix expansions: 


b(s, t, X, y) - Pe[s,t,x,y)\ = \'^f>® (s, t, x,y) - (s, t, x, y) \ 

r—0 r—0 

oo oo 

< \{P-Pe){s,t,x,y)\ + I'^p® H'''^'>{s,t,x,y) - '^Pe® H^'^\s,t,x,y)\. 


(3.1) 


r=l 


r—1 


The strategy to study the above difference, using some well known properties of the Gaussian kernels and their 
derivatives recalled in (12.81) . consists in first studying the difference of the main terms. 

We have the following Lemma. 

Lemma 1 (Difference of the first terms and their derivatives). Under (A), there exist Ci > 1, c £ (0,1] s.t. 
for all 0 < s < t, (x, y) £ (R.'^)^ and all multi-index a, jaj < 4, 


\Dfp{s,t,x,y) - D2p^{s,t,x,y)\ < 


Cl 


(t- s)l“l/2 

Proof. Let us first consider jaj = 0 and introduce some notations. Set: 


^s,<J,-yPc{t - S,y- X). 


^{s,t,y)-= a{u,y)du,®,e{s,t,y)-= ae{u,y)du. 

J s J s 


(3.2) 


Let us now identify the columns of the matrices S(s, t, y), ^^(s, t, y) with d-dimensional column vectors, i.e. for 
E(s,t,?/): 

E(s,t,y)= ( £1 I S2 I ... I )(s,t,y). 

We now rewrite p{s, t, x, y) and ps{s, t, x, y) in terms of vectors 0(s, t, y), 0e(s, t, y) £ : 

p(s,t,x,y) = /a:,y(0(s,t,y)), 0(s,t,y) = ((E^)*,... , (E‘')*)*(s, t, y), 

Pe { s , t , x , y ) =/:r.y(0e(s,t,y)), 0e(s,t,y) = ((E^)*,... , (E^)*)*(s,t,y), 

where (.)* stands for the transpose and 

fx,y '■ K 

r ^ fx.v{T) = 


exp (' ) {y--^) ] ^ 


where L := 


/ ri \ 
r2 


V rV 


(27r)d/2det(ri^'^)i/2 

and each (r®)ig|i belongs to Also, we have denoted: 

ri:<i := ( pi I r2 I ... I ), 


(3.3) 


the d X d matrix formed with the entries (ri)jg|ieach entry being viewed as a column. 
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The multidimensional Taylor expansion now gives: 


\{p-Pe){s,t,x,y)\ = |/a:,y(0(s,<,l/)) - f^^y{Qe{s,t,y))\ 


kl=2 


= XI t, i/)){(0e - 0)(s, t, y)Y 

kl=i 

/ {l-6)D''f,,yi[Q + S{e,-Q)]is,t,y))d6 , 
Jo 


(3.4) 


where for a multi-index v := (ui, • • • , Vj2) £ , we denote by \v\ := length of the multi-index, 

= nf=i and for /i e h’' := nf=i (with the convention that 0° = 1). With these notations, from 
(|3.2D . (13.31) , (13.41) we get: 

|/x.y(0(s,i,2/)) - fx,y{Qe{s,t,y))\ < c| X \D'' fx,y{Q(.S,t,y))\I\e,a,^{t “ s) 

^ IH=1 

+ ^£.<T. 7 (^-s)^ max X \D'' fx,yi[^ + - Q)]{s,t,y))\V (3.5) 

^ ’ > 1=2 ^ 

recalling that has previously been defined in dniQ. 

Since fx,y in (13.31) is a Gaussian density in the parameters x, y, we recall from Cramer and Leadbetter [CL04) 
(see eq. (2.10.3) therein), that for all T G and any multi index u, \v\ < 2: 


D“i=..,x )= 5 m n( 


92 




-rixA^) 


:-L- 


ijd 


where [-J stands for the integer part. Hence, taking from (13.51) . for all <5 G [0,1], Ti;^s{s,t,y) := [0 -I- d(0£ — 
0 )](s,t,?/) yields, thanks to the non-degeneracy conditions (see equation (12.81) '): 


\D’'fx,y{T^^s{s,t,y))\ < 


Cl 


[t — s)l'^l 


fx,y(^cTgg(^S,t,y)'j < 


Cl 


{t — s)l'^l 


Pc{t- s,y-x), 


(3.6) 


for some ci > 1, c G (0,1]. 

Thus, from (13.3|) . (13.41) . equations (13.51) and (13.61) give: 

\p{s,t,x,y) - Peis,t,X,y)\ < Cil^e,o,-iPc{t - S,y- x). 

Up to a modification of ci, c or ci, c in (12.81) we can assume that the statement of the lemma and (12.81) hold with 
the same constants ci, c. The bounds for the derivatives are established similarly using the controls of (12. 8 |) . 
This concludes the proof. □ 

Remark 2 . Observe from equation (13.41) that the previous Lemma still holds with JS.e,iy,-y replaced by Ag^CT.oo ■= 
supjgjQ |o'(t,.) — CFsit, .)|oo. The Holder norm is required to control the differences of the parametrix kernels. 

The previous lemma quantifies how close are the main parts of the expansions. To proceed we need to 
consider the difference between the one-step convolutions. Combining the estimates of Lemmas [1] and [5] below 
will yield by induction the result stated in Theorem |TJ 

^Actually only the sup norm term in could be considered for this part of the anlalysis. 
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Lemma 2 (Control of the one-step convolution). For all 0 < s < t < T, {x, y) € and for q S (d, -l-oo] ; 

\p® -pe® H'^^\s,t,x,y)\ 

< c\pc{t - s, y - a;)|2(l V + Iq=+ooAe,b,+oo].B(l, ^)(t - s)i (3 7 ) 

+IIgG(d,+oo) a{q), a{q)){t - |, 

where c\, c are as in Lemma[J\ and for q S (d, -|-oo] we set a{q) = ^(1 — The above control then yields for a 
fixed q € (d, -|-oo]; 


\p® H^^'>{s,t,x,y) - pe® HP{s,t,x,y)\ 

< 2C'^A^^^^qPcit - s,y - x){t - V B{^ + a{q), a{q))), C = ci(l V 


(3.8) 


which will be useful for the iteration (see Lemma 0). 

Proof. Let us write: 

\p® H^^1{s,t,x,y) -pe® H^^\s,t,x,y)\ < \p - Pe\ ® \H\{s,t,x,y) +Pe®\H - Hs\{s,t,x,y) =: I + II. (3.9) 
From Lemma [T] and (12.91) we readily get for all q G (d, -l-oo]: 

\p-Pe\® \H\{s,t,x,y)\ < {{IV - s,y - x)B{l,'^){t - s)i . (3.10) 

Now we will establish that for all 0 < u < t < T, (z, y) G (R"^)^ and q = -l-oo: 

fi V r(^“'>')/^')ci 

\{H - Hs){u,t,z,y)\ < As,.y,oo —--- pc{t-u,y-z). (3.11) 

{t — u) 2 

Equations (13.111) and (12.81) give that II can be handled as I which yields the result for q = -|-oo. It therefore 
remains to prove (|3.1ip . Let us write with the notations of (13.3^ : 


{H - He){u,t,z,y) := 


-Tr( {a{u,z) - a{u,y))Dlf^^y{Q{u,t,y)) ) -h {b{u, z), DJz^y{e{u,t,y))) 


^Tr( {ae{u,z) - as{u,y))Dlf^^y{Qs{u,t,y)) ) + {bs{u, z), D^f^^y{Qe{u,t,y))) 


Thus, 


{H - Hs){u,t,z,y) = 


1 


+ 


Tr^(a(u,z) - a{u,y)){Dlf^^y{Q{u,t,y)) - (©^(u, t, i/))} 

-Tr^[(ae(u,z) - ae{u,y) - {a{u,z) - a{u,y))]Dlf^^y{Qs{u,t,y)) 
{b{u,z),{DJ;,^y{Q{u,t,y)) - (©^(u, t, y))}) - {{be{u,z) - b{u, z)), DJ;,_y{Qe{u,t,y))) 


(3.12) 
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Observe now that, similarly to (13.61) one has for all i G {1, 2}: 


\Dlfz,y{Q{u,t,y))\ + - u,y - z), 

ciAe 


\Dlfz,y{Q{u,t,y)) - Dlf^^y{eeiu,t,y))\ < 


{t — uy/"^ 


ps{t -u,y- z). 


Also, 


\{ae{u,z) - ae{u,y) - {a{u,z) - a{u,y))\ < - y^, 

|6e(m, z) - b{u, z)\ < cAs^b,oo- 

Thus, provided that ci, c have been chosen large and small enough respectively in Lemma [U the dehnition in 
(cni) gives: 


\{H - H^){u,t,z,y)\ < - (^_^)i- 7/2 Pc{t-u,y-z). 

This establishes (13.1111 for q = +oo. For q G (d, +oo) we have to use Holder’s inequality in the time-space 
convolution involving the difference of the drifts (last term in (13.121) '). Set: 

D{s,t,x,y) := du Pe{s,u,x,z){{be{u,z)-b{u,z)),Dzfz^y{<ds{u,t,y)))dz. 

J s Jm'^ 


1/9 


Denoting by q the conjugate of g, i.e. q,q > l,q ^ + q ^ = 1, we get from (12.811 and for g > d that: 

\D{s,t,x,y)\ < “ ^e('“>-)llL-!(R<i){ J J,Pc{u - s, z - x)pc(,t - u,y - z)ydz''^ 

d~^^s,b,qPcit-s,y-x) 


< ci 


du 


( u — s )'^^^ — 11)2 + 2(1 


Now, the constraint d < q < -l-oo precisely gives that 1 < g < d/(d — 1) ^ -I- |(1 — |) < 1 so that the last 

integral is well defined. We therefore derive: 

\D{s,t,x,y)\ < cfit - s,y-x)B{l - ^(1 “ t), ^ - ^(1 - |)). 


In the case d < q < -foo, recalling that a;(g) = i(l — |), we eventually get : 


Pe (g) |7d - Hy{sy,x,y) < c\pc{t - s,y- x){As,b,q(t - + a{q),a{q)) 

+2Ae.<,,^(i V - sy/^Biy 7 / 2 )}. 


(3.13) 


The statement now follows in whole generality from (13.9p , (j3.10|l , equations (13.111) , p. 8 ll for g = oo and (I3.13p 
for d < g < - 1 - 00 . □ 

The following Lemma associated with Lemmas [1] and [2] allows to complete the proof of Theorem [1] 
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Lemma 3 (Difference of the iterated kernels). For all 0 < s < t < T, {x,y) G and for all q G (d,+oo], 

r e N; 




Pe 0 H^^'>){s,t,x,y)\ < (r + 


(7-+1 [r(lAa(g))]’~ 
T{1 + r{l A a{q))) 


Pc{t- s,y 




(3.14) 


where c is as in Lemma [T] and C as in equation (13.81) of Lemma [2] 

Proof. Observe that Lemmas [T] and O respectively give (13.141) for r = 0 and r = 1. Let us assume that it holds 
for a given r G N* and let us prove it for r + 1. 

Let us denote for all r > 1, Pris, t, x, y) := |(p 0 — Pe 0 F[e^'^){s, t, x, y)\. Write 

Pr+iis,t,x,y) < |[p0 -Pe 0 0 H{s,t,x,y) \ + \pe 0 0 (id - Hs){s,t,x,y)\ 

<Pr® \H\{s,t,x,y) + \pe 0 hP\ 0 \{H - He)\{s,t,x,y). 


Recall now that under (A), the terms \F[\{s,t,x,y) and \pe 0 iii^^| satisfy respectively and uniformly in e the 
controls of equations (12.91) , ( 12 . 101 ) . The result then follows from the proof of Lemma [5] (see equation (13.111) for 
q = oo and (13.131) for q G (d, +oo)) and the induction hypothesis. □ 

Theorem [I] now simply follows from the controls of Lemma [31 the parametrix expansions (12.61) and (12.71) of 
the densities p, Pe and the asymptotics of the gamma function. 


3.2. Stability for Markov Chains. 

In this Section we prove Theorem |3J The strategy is rather similar to the one of Section 13.11 thanks to the 
series representation of the densities of the chains given in Proposition 

Recall first from Section ISJ that we have the following representations for the density p^,p^ of the Markov 
chains Y, in (11.31) . For all 0 < ti < tj < T, {x, y) G (R"^)^: 


j-i 

p’^{ti,tj,x,y) = ^p^ 0?i H^'^"’\ti,tj,x,y), 

r—0 

j-i 

Pe{U,tj,x,y) = ^p^ H^’^'"Hu,tj,x,y). 

r—0 


3.2.1. Comparison of the frozen densities 

The first key point for the analysis with Markov chains is the following Lemma. 

Lemma 4 (Controls and Comparison of the densities and their derivatives). There exist c, ci s.t. for all 
0 < ti < tj < T, (x, y) G (R'^)^ and for all multi-index a, |a| < 4; 




|D“p (t*,tj,x,y) - D“p^(ti,tj,x,j/)| < 




{t,-U)\<^\l^ 


^C,C1 (fj 

{ij 


-ti,y-x), 

-ti,y-x), 


where 

- Under (Ig)-' 

i’c,ciitj -ti,y- x) := ciPcftj -ti,y - x), 
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- Under (Ip^m): 

Proof. Note first that under (Ig) the statement has already been proved in Lemma [T] We thus assume that 
(Ip,m) holds. Introduce first the random vectors with zero mean: 






1 


'^(7e{ti,y)Vh£,l+l. 

l=k 


Denoting by qj-k,qj-k,e their respective densities, one has: 

1 


Dfp (tk,tj.x,y) = 
^xPei^kUjj^^y) — 




(3.15) 


From the Edgeworth expansion of Theorem 19.3 in Bhattacharya and Rao [BR76) . for qj-k,qj-k,e^ one readily 
derives under (A), for |a| = 0 that there exists ci s.t. for all 0 < tk < tj < T, {x,y) S 


p'"{tk,tj,x,y) +p'f{tk,tj,x,y) < 


{tj 


Cl 1 



\x-y\ \ 


(3.16) 


for all integer m < M — d, where we recall that M stands for the initial decay of the density /c of the innovations 

bounded by Qm (see equation (ini))- _ _ 

We can as well derive similarly to the proof of Theorem 19.3 in [BR,76] . see also Lemma 3.7 in [KMOO] . that 
for all a, |a| < 4: 


DxP (tk,tj,x,y)\ + \D^p^{tk,tj,x,y)\ < 


Cl 


{tj - 4)(''+l“l)/2 , \x-y\ y 

V / 


(3.17) 


for all m < M — d — 4. Note indeed that differentiating in Df the density and the terms of the Edgeworth 
expansion corresponds to a multiplication of the Fourier transforms involved by C standing for the Fourier 
variable. Hence, from our smoothness assumptions in (Ipm), after obvious modifications, the estimates of 
Theorem 9.11 and Lemma 14.3 from [BR,76] apply for these derivatives. With these bounds, one then simply 
has to copy the proof of Theorem 19.3. Roughly speaking, taking derivatives deteriorates the concentration of the 
initial control in (13.161) up to the derivation order. On the other hand, the bound in (13.161) is itself deteriorated 
w.r.t. the initial concentration condition in dni. The key point is that the techniques of Theorem 19.3 in [BR76) 
actually provide concentration bounds for inhomogeneous sums of random variables with concentration as in 
dm in terms of the moments of the innovations. To explain the bound in (13.161) let us observe that the m**' 
moment of ^ is finite for m < M — d. 

Equations (13.161) and (13.171) give the first part of the lemma. Still from the proof of Theorem 19.3 in |BR76] . 
one gets, under (A), that there exists C > 0 s.t. for all multi-indexes a, |q:| < 4,,3, |,9| < m < M — d — 5 for 
all j > k: 

\C\{\Dfq,.kiC)\ + \Dfq,.kAO\yC < C, (3.18) 

where qj-k{C),qj-k,eiC) stand for the respective characteristic functions of the random variables j, at 

point Q. 
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a\ < 


To investigate the quantity y) — D^p^{tk, tj,x, y)\ thanks to (13.151) define now for all a, ^ _ 

4,/3, \/3\ < m < M — d — 5: 

Vz e Qj-k,e{z) := {qj-k{z) - qj-k,s{z )), 

VC e Qj-kAC) ■■= (C“ {qj-kiO - qj-kAC)}) ■ (3-19) 

Let us now estimate the difference between the characteristic functions. From the Leibniz formula, we are 
led to investigate for all multi-indexes /3,d, \^\ < |/3|, |q;| < |a| quantities of the form: 

= exp[*C • ZIJ - {Zl^^f exp[*C • Zl^^_ . 

Assume first that j > k + \. In that case, set now Zf, - ^ := j 2 j ~ j i- Denoting 

similarly := ■= perturbed process, we get: 

{z^^)-^C{Dlq,.k[0 - Diq,_kAO) = 


C“< E 


(^fcj.i + ^h,2f exp[< • ZA^i] exp[< • ZA^2\ “ ^ + Zlf^f exp[< • Zlf}\ exp[< • 


I Z.|i|<|^| 

^ C^E expK • Zl^^h] E [(Zliif-^ exp[*C • | = 

E ® ®^P[< • ® exp[*C • Zl, , 

I /,|i|<|^| 

(^Si exp[< • Zl’j'h] [e[{ZI^A^~^ ®^pK ■ ^1^4 exp[zC • Zl’jy 


E 


where in the above expression we considered the binomial expansion for multi-indexes denoting by := 
with the corresponding definitions for factorials (see the proof of Lemma [T]). Introduce now, for a multi-index 
I, |Z| S |0, |,5|], the functions: 


CE[iZA ,f-^ exp[zC • 4",,- 2 ]], ^/“’'(C) := C'*E[(%g)' exp[*C • 


and 


SuiO ■■= [E[izy^^A'^^pl<-^h4 -4^zif}yeMK-zlfl]_ 

:= [E[izy^^A^-AM<-zA,,]] -E[{zyf}f-AMK-~zli%_ 

Thus, we can rewrite from the previous computations: 

(yA-H^{Dlq,.k{0 - Dlq^-kAO) = E {(£:i,z4-f’'^-')(C) + (^2. 

i.|Z|<|^| 


(3.20) 

















































18 


STABILITY OF DENSITIES FOR PERTURBED DIFFUSIONS AND MARKOV CHAINS 


Recall from (I3.18|) that we already have integrability for the contributions and d> 2 ’^(C)- Let us thus 

proceed with the control of We only give details for the contribution £ 2,^-1 can be 

handled similarly. We also consider |Z| > 2, since the cases |?| < 2 can be handled more directly. Write: 






< C 


{ nzl ,,. 


^kj,l 


\i\zi. 


fej'T I 




\z> 


vAs) I 
fcj ',1 I 


I^ 11 (C)I < 


exp« • 1 ) - exp« ■ 1 

dSi'‘iciizL,, - 


')| + E||Z: 


1 }^ 


Apply now Holder’s inequality with Pi = \l\,qi = |Z|/(|^| —1) for the first term and p 2 = (|^| + l)/|^|j 92 = |/| + lfor 
the second one so that all the contribution appear with the same power (in order to equilibrate the constraints 
concerning the intregrability conditions). One gets: 


|C|E[|Z^;g|l'l+i]l'l/(l'l+i)]E[|^y^. ^ _ ^yjO||J|+i]i/(|/|+i)|^ (3 21) 

The point is now to prove, since we have assumed m < M — d — 5 m + 1 < M — d — 4, that there exists c 
s.t. for all r < TO + 1 , 


n\zi,,,-zif£!ry/^<cA, 


£,< 7,75 


mi 


fcj'.i I 


ll/r 


■E[|Z 


fej'.l I 


ll/r 


< C. 


(3.22) 


Let us establish the point for the difference, the other bounds can be derived similarly. Define for all i G 
|fc, j]. Mi := VhJ2l-Jk{^ ~ o'e){tr, y)^r+i- The process {Mi)i^ik,j} is a square integrable martingale (in discrete 
time, w.r.t. Ti := < i), S-field generated by the innovation up to the current time). Its quadratic 

variation writes [M]i = K'^ ~ o’E)(^r, 2 /)Pl^r+iP and the Burkholder-Davies-Gundy inequalities, see e.g. 

Shiryaev [Shi96) . give for all r < M — d — 4: 


i-i 


E[ sup \M^l] < CrE[[M]j/^] = Cr/i’'/^E[(V |(cr - as){ti,y)y\Ci+iyY^'^]- 


i—k 


Hr = 2 one readily gets: 

nzL .,- zim ]< 


j-i 




E[ sup < 


{tj tk) iG[fcj] 
Let us thus assume r > 2 and derive from (13.2311 


{tj tfe) 


7' 


i—k 


nzi,,i - zi:f}n < %y,2 n sup m 

{tj-tk)' jG[fcj1 


< 




J -1 


1-1 


[tj - tkY'^ 


e[(E - ^E)(iDj/)ri?.+ir)(E 


(3.23) 


applying Holder’s inequality for the counting measure with p = r/2,(; = rjir — 2) for the last inequality. This 
finally gives: 




1-1 


E[l%i - zlfjY ] < 
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Since we have assumed r < m + 1 < M — d — A, this gives the first control in (13.221) . The other one readily 
follows replacing a — hy u oi 

From equations (I3.2ip . (I3.22p and similar controls for 82 p-iiO we finally derive: 

l^fo(C)l + |f2j-l(C)l <t^lAe,..^(l + |C|). 

As a result we have from (13.1911 and (j3.20l) : 

E (i^?’''”'(c)i + i'i^2’'(c)i)(i + ici) 

l ^.l/J|<l/3L iTI<l/3| 

a = a - {P - P). 


We finally derive from (13.191) and (13.181) 
as well for a multi-index a, |d| = 5): 


(which thanks to the smoothness assumption on Qm in (Ip,m) holds 




< 





(3.24) 


From (13.1511 this concludes the proof for j > k 1- li j = k + \ the previous arguments can be simplified and 
lead to the same results. □ 


3.2.2. Comparison of the parametrix kernels 

This step is crucial and actually the key to the result for the Markov chains. We focus for simplicity on the 
case q = -boo, for which pointwise controls for the differences between the drift coefficients are available, and 
which already emphasizes all the difficulties. The case q G {d, + 00 ) for the drifts could be handled as in Lemma 
[21 using similar Holder inequalities. 

We actually have the following Lemma. 


Lemma 5 (Control of the One-Step Convolution for the Chain.). There exists ci,c s.t. for all 0 < tk < tj < 
T,iz,y)G(R^)^: 


\{H^ 


H^){tk,t,,z,y)\ < 


tk,z- y), 


with 

- -tk,z - y) = ipc,ci{tj -tk,z-y) under (Ig). 

- $c,ci [tj - tk,z - y) = ifc,ci (tj - tk, z - y) (1 + ’ under (Ip,m), 

where ipc,ci is defined according to the assumptions on the innovations in Lemma^ 


Proof The case k = j + 1 involves directly differences of densities and could be treated more directly than the 
case fc > j + 1. We thus focus on the latter. Introduce for k G |0, N], {x, w) G (M'^)^ the one step transitions: 


w, 


T^{tk,x,w) := b(tk,x)h + h^/^a{tk,x)w, T^(tk,x,w) := be{tk,x)h + h^/‘^ae{tk,x) 
TQ{tk,x,w) := h^/'^a{tk,x)w, T^Q(tk,x,w) := h^^‘^ae(tk, x)w. 


(3.25) 
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From the definition of recalling that /j stands for the density of the innovation, the difference of the 

kernels writes: 


{H^ - f dwf^iw) z + z,w),y) - z + Tl^{tk,y,w),y)\ 

jRd '' > 


- {pe (4+1 ,4, z + (tfe, z, w), J/) - {tk+1 ,tj,Z + (tfe , j/, it) , y) I 


(3.26) 


Let us now perform a Taylor expansion at order 2 with integral rest. To this end, let us first introduce for 
A G [0,1] the mappings: 




V’t'- 


(ri,T2) ^ Tr(i5,V(4+i4„^ + ATi,y)[T2T2*]), 

X K 


(3.27) 


(Ti,T 2) ^ Tr(i?2p^^(4+i,t„z + ATi,y)[T2T|]), 

where T 2 is viewed as a column vector and T| denotes its transpose. Recalling as well that ^ is centered we get: 

:= z,y) 

+h~^ f dwf^{w) f dX{l — A) 

Jo 

z, w), T^{tk,z, wj) - (px{TQ{tk,y, w), T^itk, y, w))| 

-{pA.e(^^(4, z, w), T^itk, z,w))- w), T^^^itk, y, w))| 


=: (AiiL'‘'^ + A2iL'‘’^)(4,t„z,y), 


(3.28) 


where for i G {1,2}, is associated with the terms of order i. The idea is now to make appear 

explicitly. The term AiiL^’^ is the easiest to handle. We can indeed readily write: 


AiH^’^tk,tj,z,y) 

(^D^p^{tk+i,tj,z,y), [b{tk,z) - be{tk,z)]j - (^{D^p^ - D^p^){tk+i,tj, z,y),be{tk, z)'j 

From Assumption (A3), equation (11.51) and Lemma 0] we derive for q = +oo: 

CA, 


AiiL'‘’"(tfc,t„z,y)| < 


i£,7,00 


(4-4)42 


i^c,ci{tj -tk,y- z). 


(3.29) 
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The term is trickier to handle. Define to this end: 


Aip’l’‘'{tk,z,y,w) := z,w),T^{tk, z,w)) - ipl{Tg{tk,y,w),T^{tk,y,w))'^ 


Let us then decompose: 


A(p’l'%tk,z,y,w) := z,w),T^(tk, z,w)) - ip’KT^itk, z,w),TQ{tk,y,w))^ 

-{A,eiT^itk,z,w),T^{tk,z,w)) - ip'l^^{T^{tk,z,w),Tl^^^{tk,y,w))^ 
[‘P\{T''{tk,z, w), T^{tk,y, w)) - ip'l{T^{tk, y, w),T^{tk, y, w))| 
-{v>X,eiT^iik,y, w), T^^^itk, y, w)) - <^A.e(^o!e(^fe: w), T^^^(tk,y, w))| 


+ 


=: + ^ 2 ‘Px'^){tk, z, y, w) 


h , e \ 


(3.30) 


and write from (j3.28|) : 


A 2 H’^''^{tk,tj,z,y) = h ^ f dwf^{w) ( (iA(l - A)(Ai(/?^’® + A2V3 a ^)(4 , 2 :, y, w) 

Jo 

=: {A 2 iH^’^ + A 22 H^'^){tk, tj, z, y), 


(3.31) 


for the associated contributions in A 2 H^'^. Again, we have to consider these two terms separately. 
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Term We first write from (13.301) : 


z, w), T^(tk, z, w)) - z, w), T^itk,y, w))| - 

{v>xiT^itk,z, w),T^{tk, z, w)) - ipx(T^itk, z, w), T^^^{tk,y, rt;))| 


{vx{T^{tk,z, w), T^{tk, z, w)) - z, w), T^^^{tk,y, w;))| 


-[y^'liT^(tk,z,w),T^{tk,z,w)) - ipl(T^{tk,z,w),T^^^(tk,y,w))'^ 
[‘Px,6iT^itk,z,w),T^{tk,z,w)) - ip'l ,^{T^{tk,z,w),T^,^{tk,y,w))'^ 
-[v>xiT^(tk,z,w),T^{tk,z,w)) - ip’]i(T^{tk,z,w),Tl^^^(tk,y,w))'^ 


3 

= : '^Au(p'l'‘'{tk,z,y,w). 
2=1 


We now state some useful controls for the analysis. Namely, setting: 

D{tk, z, y, w) := r^(tfe, z, w)T^{tk, z, w)* - T^(tk,y, w)T^{tk, y, w)*, 
De{tk,z,y,w) := T^(tk, z,w)T^(tk, z,w)* - Tg^^(tk,y,w)TQ^^{tk,y,w)*, 

we have from (A3) and equation p.5l) for q = +oo : 

{\D\ + \De\){tk,z,y,w) < c(/i^ ++/i(l A |z - y|)'>'|w|^), 

\D - Di;\{tk,z,y,w) < + h^^'^\w\ + h{lA\z-y\y\w\'^). 

From the definition of (p^ in (13.271) . equation (13.321) . the control (13.331) and Lemma 01 we get: 


\‘^ii>fx^\{tk,z,y,w) 

< - ^ ^ 3 ,, 1^1 ^ ^ 

[tj tk) 


(3.32) 


(3.33) 


(3.34) 
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We would similarly get from Lemma S] and (13.331) : 


< rA 
— *^^£, 7,00 


4’c,ci (7 -tk,y - jz + XT^jtk, Z, w))) 
(tj — tk) 


+ h^/^\w\ + h{lA\z-y\y\w\^), 

\^i2V>\^\itk,z,y,w) 


< 


V^c,ci (7 


< cA, 


V^c,ci (ij 


— ‘-^£, 7,00 


tk,y- {z + dXT’^{tk,z,w) + (1 - 9)XT^{tk,z,w))) 
{tj - tkf/^ 

x|(r^ - T^){tk.,z,w)\\De\{tk,z,y,w) 
tk,y - (z + 9\T^{tk,z,w) + (1 - 0)XT^{tk,z,w))) 
{t^-tkf/^ 

x{h^ + h^/^\w\ + h{l A\z- y\y\w\^){h + h^/^\w\), 


(3.35) 


for some 9 S (0,1), using as well (13.251) and (11.51) for the last inequality. The point is now to get rid of the 
transitions appearing in the function V'c.ci. We separate here the two assumptions at hand. 


Under (Ig), it suffices to remark that by the convexity inequality \z — y — Q\‘^ > i|z —yp — |0p, for all 0 G 


V'c.ci(7 - tk,y - z - Q) < Cl 


rd /2 


(27r(tj 


exp 


c \z-y\^ 

4 tj - tk 


exp 77 


| 0 | = 


2 tj tk 


Now, if 0 is one of the above transitions or linear combination of transitions, we get from (13.251) : 


V’c.ci {tj -tk,y - z - Q) < Cl 


{cl2)d/2 


(27r(7 - tk))d/'2 


exp - 


c \z-y\^ 

4 tj - tk 




(3.36) 


up to a modification of ci observing that h/{tj — tk) < 1 and with K 2 as in (Ai). Since c can be chosen small 
enough in the previous controls, up to deteriorating the concentration properties in Lemma SI the last term 
can be integrated by the standard Gaussian density appearing in (I3.31|) . We thus derive, from (13.361) . (I3.34|) . 
(j3.35|) and the definition in (j3.32l) , up to modifications of c, ci: 


\^i^x‘'\{tk,z,y,w) < 


Ae,^,oohCT/ic,ci {tj -tk,z-y) exp(c|r(;|^) < 1 + 




— yl'^l'7/Tp 


{tj - 4)^/^ tj - tk 


which plugged into (13.311) yields up to modifications of c, c, Ci: 

Ae,.^,oo(l V T(^“'>'^/^)-0 c.ci(7 -tk,z- y) 


|A 2 li?^’ 7 tfe, 7 , 2 :, 1 /)| < c- 


(7 - tk)^-'^/^ 


(3.37) 


- Under (Irm), we only detail the computations for the off diagonal regime \z — y\> c{tj — tk)^!"^ which is the 
most delicate to handle. In this case, we have to discuss according to the position of w w.r.t. y — z. With the 
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notations of (A2), introduce := {w G : {A/i}^/^|w| < 1^: — ?/|/2}. liw&V, then, still from (13.341) . (I3.35p . 


< — —{h"^ + h^^'^\w\ + h{l A\z- y\y\w\'^), 

[tj tk) 

\^i 2 V’x^\{tk,z,y,w) 

< cAe,^,oo + h^/‘^\w\ + h{l A\z- y\y\w\‘^){h + h^/'^\w\). 

On the other hand, when w ^ V we use to make the off-diagonal bound of '<pc,ci{tj — tk,y — z) appear. 
Namely, we can write: 


/cH 


< 


< 


1 „ 1 1 
^(1 -f |w|)*f “ ^(1+ \^-y\ )M-(d+i) (l-p|w|)^+4 
1 1 
""(I + + 1^1)^+^’ 


(3.38) 


where the last splitting is performed in order to integrate the contribution in |?iip coming from the upper bound 
for \Ai 2 Lp^^\ in (13.3511 . Plugging the above controls in (I3.31|) yields: 


AH^{^{tk,tj,z,y) \ < 


Ag^^-jOO^c,ci ik: Z y) 


(3.39) 


We emphasize that in the case of innovations with polynomial decays, the control on the difference of the kernels 
again induces a loss of concentration of order 7 in order to equilibrate the time singularity. 


Term A 22 H^’^- This term can be handled with the same arguments as A 2 iH^’^■ For the sake of completeness 
we anyhow specify how the different contributions appear. Namely, with the notations of (13.301) and (13.311) : 


A2(Px''itk,z,y,w) = 


dfj.^(^DT,ipx{To{tk,y,w) + fi(T’^{tk, z,w) - TQ{tk,y,w)),TQ{tk,y,w)),T^{tk,z,w) - T^{tk,y,w) 
Dr^^xA'^oA^k.y.w) + ^i{T^{tk,z,w) - ro'(^(tfc,y,u;)),To^^(tfe,y,u;)),Tg^(tfc,z,u;) - To'(^(tfe,y,u;)^| 
= {/ d^i^(pT^Lp\{TQ{tk,y,w) + y,(T^{tk,z,w) - Tg{tk,y,w)),Tg{tk,y,w)), 
[{T^{tk,z,w) - TQ{tk,y,w)) - {T^(tk,z,w) - TQ^^{tk,y,w))]'j'^ 


dfj. 


\DT^V\,eiTo,s(dk,y,w) + y.{T^{tk,z,w) - T^^^{tk,y,w)),T^^^(tk,y,w)) 
-DTi‘pl{Tg{tk,y,w) + fi(T^{tk,z,w) - TQ(tk,y,w)),TQ{tk,y,w)) ,T^{tk,z,w) - TQ,^(tk,y,w)^'^ 

= ■■ (A21V?a ^ + ^22V’x^)itk, Z, y, w). 
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In A 2 iip^’^ we have sensitivities of order 3 for the density, giving time singularities in {tj — tk) which are 
again equilibrated by the the multiplicative factor: 


\Toitk,y,w)[T^{tk,y,w)]*\ 
x\{T'^(tk,z,w) - T^{tu,y,w)) - {T^{tk,z,w) - T^,^{tk,y,w))\ 

< c{h^ + A\z- yD'^lwl), 

where the last inequality is obtained similarly to (13.331) using as well (11.51) . The same kind of controls can be 
established for A 22 <P\’’^- Anyhow, the analysis of this term leads to investigate the difference of third order 
derivatives, which finally yields contributions involving derivatives of order four. This is what induces the final 
concentration loss under (Ip^m), i-e. we need to integrate a term in (see also equation (I3.38P in which we 
performed the splitting of /c on the off-diagonal region to integrate a contribution in |wp). 

We can thus claim that 


\AH^^fftk,tj,z,y)\ < 


^£,7,00^c,Cl y) 

(7 - ■ 


Plugging the above control and (j3.39|) (or (13.371) under (Iq)) into (j3.31|) we derive: 


\AH!^’’'{tk,tj,z,y)\ < 


^c,ci (^^' y) 


which together with (13.291) and the decomposition (13.281) completes the proof. □ 

From Lemmas [4] and [5] the proof of Theorem [2] is achieved, under (Iq), following the steps of Lemmas [2] and 
[21 using the Holder inequalities for the differences of the drift terms for q G (d, +oo). 

The point is that we want to justify the following inequality under (Ip^m) and q = +oo: 


< (r + 


{(1 [r(§)]’ 


|(p^ -pg ®h H^'^''^){U,tj,x,y)\ (3.40) 

y-x 


r(i + 


ri) 


(7 - 


Qm—{ tZ+5+7) 


(7 - uff/^/c 


{tj -ti)^ . 


The only delicate point, w.r.t. the analysis performed for diffusions, consists in controlling the convolutions of the 
densities with polynomial decay. To this end, we can adapt a technique used by Kolokoltsov [KolOO] to investigate 
convolutions of “stable like” densities. Set to := M — (d + 5 + 7 ) and denote for all 0 < i < j < A, a: G by 

qm{tj - U, x) := (^^._‘).)d/ 2 QM-(d+ 5 + 7 ) ( fa_tf)i/ 2 /e ) the density with polynomial decay appearing in LemmasH 
andjSJ Let us consider for fixed i < k < j, {x, y) G (R'^)^ the convolution: 


ll^{ti,tj,x,y) := / dzqm{tk - ti,z - x)qm{tj - tk,y - z). (3-41) 

Jr‘‘ 

- If |a: —y| < c{tj (diagonal regime for the parabolic scaling), it is easily seen that one of the two densities 

in the integral (13.411) is homogeneous to qm,{tj —ti,y — x). Namely, if {tk — U) > {tj — ti)l2, qm{tk — ti,z — x) < 

^ A cqm{tj ti,y x). Thus, 

4 {u, tj,x, y) < cqm{tj -u,y-x) f dzq^tj -tk,y-z) = cq^{tj -U,y- x). 

dR<J 

If {tk — ti) < {tj — ti)/2, the same operation can be performed taking qm{tj — tk,y — z) out of the integral, 
observing again that in that case qm{tj — tk,y — z) < cqm{tj — ti,y — x). 
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- li\x — y\> c(tj — (off-diagonal regime), we introduce Ai := {z G : |a: — 2;| > ^\x — y\}, A2 '■= {z G : 

\z — y\ > \\x — y\}- Every z G belongs at least to one of the {^i}iG{i, 2 }- Let us assume w.l.o.g. that z G A 2 . 
Then \z — y\ > > ^{tj — so that the density qm{tj — tk,y — z) is itself in the off-diagonal 

regime. Write: 


dzq^itk - U,z - x)qm{tj - tk,y - z) < / dzqmitk - U,z - x) 


Cm(tk - 

\z - y|™ 


< 


\x - y\^ 


dzq^{tk -U,z-x) < cqrnitj -ti,y- x), 


recalling that, under (Ip,m)) "m- > d for the last but one inequality. The same operation could be performed on 

Ai- 

We have thus established that, there exist c > 1 s.t. for all 0 < i < fc < j, {x, y) G : 


{ti ,tj,x,y) < cqm {tj -U,y-x). 

From the controls of Lemma [5] and following the strategy of Lemma [31 we will be led to consider convolutions 
of the previous type involving T functions. The above strategy thus yields (13.401) by induction. 
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